Discontinuity in the specific heat of a weakly interacting Bose gas 
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We produce the discontinuity in the specific heat of a homogeneous, dilute, and weakly interacting 
Bose gas in a short-wavelength range with a simple statistical method. The magnitude of the 
discontinuity at the phase transition temperature is obtained as a function of the density and 
scattering length of the Bose particles. 
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^He becomes superfluid below 2.18K, so called A point, 
at low atmospheric pressure. The specific heat of liquid 
He is infinitesimal around OK, and rises with until 
it reaches a the A point that shows a dramatic divergence, 
and then decreases tending asymptotically to a constant 
classical value. The transition is known as the strongest 
degeneracy effect of a boson system, but it is still not fully 
understood. Many physicists have tried to reproduce the 
discontinuity of the specific since the discovery of the 
transition, but even until now very little researches has 
done it successfully. 

Kikuchi et al. applied a partition function proposed 
by Feynman to an Ising model-like 2D cubic lattice and 
considering nearest-neighbors only, they calculated the 
specific heat at Tc numerically [1]. Although they pro- 
duced the discontinuity of the specific heat, their method 
strongly depends on various atomic bond models and ex- 
perimental variables. Therefore, this semi-empirical and 
numerical method could not give any formula between 
the A transition and the interaction strength. 

The objective of this paper is to reproduce the dis- 
continuity of the specific heat of a homogeneous, dilute, 
and weakly interacting (HDWI) Bose system analytically. 
Through the application of a simple statistical method, 
the relation between the A transition and the interact- 
ing strength will be obtained. The weakly interacting 
Bose model may prove useful to understand the strongly 
interacting liquid *He. We will write this argument in 
D-dimensions {2 < D < 3) for better applicability. 

The dilute is expressed by a dimensionless gas parame- 
ter 7 = n{l)a, where n(l) is the one-dimensional number 
density and a is the s-wave scattering length, fn a homo- 
geneous system, the density in D-dimensions is expressed 
as n{D) — n{l)^ . In a repulsive dilute gas, 7 is positive 
and much less than 1. The 7 of liquid ^He is known as 
on the order of one-half, and therefore the HDWI model 
is different from the real Bose system. 

For a momentum- independent interaction g{D), the 
mean field contribution to the self-energy is n{D)g{D). 
Then, the dispersion relation of the weakly interacting 



Bose gas can be written as 



^en + n{D)g{D), 



(1) 



where Eq — p^/2m. The g{D) is the positive coupling 
constant in D-dimensions, where 2 < 13 < 3. With 
the hard-sphere interaction, it is well-known in 3D as 
5(3) — Airh^a/m. The dispersion relation corresponds to 
a short-wavelength range in the Bogoliubov energy spec- 
trum 0. The exact form of the interaction is 2n{D)g{D) 
instead oin{D)g{D) because the two mean field contribu- 
tions to the self-energy from Hartree and Fock are equal. 
However, we will use Eq. ([T]) because the factor of 2 is 
not a key concept here. 

The grand partition function of the momentum- 
independent potential is given from Eq. ([1]) as 



Q(z,n,T) = Hr 
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where /3 = l/fc^T. The z is the fugacity given by z 



(2) 



and 1 below Tc- The fj, is the chemical potential, which 
is below the and negative above the Tc- The Zg is 
the effective fugacity given by 
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Note that < Zg < 1. As T -> T^, then Ze ^ e 
rjc- Note that < 77c < 1, and 1 for the non- interacting 
Bose system. The effective chemical potential is written 
in the same way : /Xg = — ng. It is —ng instead of zero 
below Tg. Even for a momentum- independent potential, 
it shifts the transition temperature. 

The equation of the state of the interacting system in 
D-dimensions is written from Eq. ([2]) as Q 



(4) 
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where A = ^j2T:h?- IrnksT is the thermal wavelength, and 
9s{z) = TldLi -zV^ is the Bose function. The only differ- 
ence of the formula from the non-interacting model is the 
effective fugacity Ze , and the information of the coupling 
constant g{D) is in it. 



2 






3 r 


CD 


2.5 - 


o 


2 - 


CD 
CU 


1.5 - 


o 


1 - 


Specf 


0.5 - 



Y=0.02 . 




7=0.01 / 




y=0.005 / 


A 


'y=0(ideal) ,f"''^- 




1 1 1 1 1 



0.2 0.4 0.6 0.8 

Temperature, T/Tg 



0.2 0.4 0.6 0.8 1 

Temperature, T/Tq 



FIG. 1: Effective fugacity of the weakly interacting Bose gas 
in three-dimensions. 



FIG. 2: Specific heat of the weakly interacting Bose gas in 
three- dimensions . 



Let To = 2Trh'^n{2)/mkB = 1 since n{D)i = n{2). It 
is about 5.92K for liquid "'He. In this way from Eq. (j4|) 
the relation of the transition temperature between ideal 
and interacting system is written as 



T^iD)^gnil)^T,iD)^gM = l. 



(5) 



The right-hand side is actually Tq. In this way the Tc(3) 
of the interacting systems is obtained self-consistently by 
numerical method as 



(6) 



The relation between Tc(3) and T°(3) is also written as 
the power series of the 7. The leading order in 3D is 



r,(3)~T0(3)(l + C7). 



(7) 



We may use the value obtained from the numerical sim- 
ulation as c ~ 1.3 [3, 4]. 

We obtain the internal energy in D-dimensions from 
C/(T) = -(a/a/3). In Q(z,n,T). 
At T < : 



U-(T) D D,, 



(8) 



and at T > T. 



U+{T) _ 

NkB ' 



(9) 



The internal energy is continous at Tc as 
limT^T, AC/(T) = 

The specific heat at constant volume is obtained as 
C,„ = {d/dT)yU{T). Therefore, we have C„ below and 
above the transition temperature in 3D. 
AtT <Tr : 



NkB 
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T-^g'M+i-^T-^g.M. 
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and at T > T. 



C+ _ 15 .g|(^e) _ 9 g|(ge) 

NkB ~ 4 g^{ze) Igdze)' 



(11) 



We used Eq. (O and the relation ^n(3)g(3) = 27/T and 
d-q/dT = 2777/T2. The discontinuity of the C„ at is 
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The non-interacting limit goes to the well-known results 
as lim^^o AC„ = 0, since 51 (1) goes to infinity. 

The effective fugacity Ze{T) = ze~^ is shown in Fig. 
1. It downshifts from 1 as the interaction is increasing. 
The Ci, at _D = 3 is shown in Fig. 2, too. 

We considered the short-wavelength range only in the 
dispersion relation, and applied a simple mean field based 
statistical mechanics to the homogeneous, dilute, and 
weakly interacting Bose system. In spite of the simplic- 
ity of the system, it successfully produced the relation 
between the discontinuity of the specific heat and the 
interaction strength of the model clearly. 
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